1 Loading of data

1.1 Data set

First, we load, filter, and merge the data sets.

How does the data set looks like

1.2 Set tresholds

Applied thresholds are indicated by grey horizontal line.

1.2.1 Area

1.2.2 Mean Intensity

1.2.3 Intensity cytoplasm ring 10

1.2.4 Intensity cytoplasm ring 9

1.2.5 Mean intensity cytoplasm ring 10

1.2.6 Mean intensity cytoplasm ring 9

1.3 Counts per sample

Cell counts per cell line:

#data <- read.csv("results/tables/data_GLUT3.csv")
table(data$Metadata_SampleID)
## 
## i1JF-R1-018 iG3G-R1-039 i1E4-R1-003 iO3H-R1-005 i82A-R1-002 iJ2C-R1-015 
##        4440        4638        4926        2816        3022        4829 
## iM89-R1-005 iC99-R1-007 iR66-R1-007 iAY6-R1-003 iPX7-R1-001 i88H-R1-002 
##        2755        3953        4291        4226        4539        4305

Mean cell count:

mean(table(data$Metadata_SampleID))
## [1] 4061.667

2 Visualize mitochondrial parameters

Various GLUT parameters are visualized for each patient-derived cell line as well as for the disease state Mean Ctrl levels are indicated by grey horizontal line.

2.1 Cytoplasmic area

2.1.1 each sample

2.1.2 disease-state

2.2 Mean intensity

2.2.1 each sample

2.2.2 disease-state

2.3 Intensity cytoplasm ring 10

2.3.1 each sample

2.3.2 disease-state

2.4 Intensity cytoplasm ring 9

2.4.1 each sample

2.4.2 disease-state

2.5 Mean intensity cytoplasm ring 10

2.5.1 each sample

2.5.2 disease-state

2.6 Mean intensity cytoplasm ring 9

2.6.1 each sample

2.6.2 disease-state

3 Statistical testing using linear mixed effects models

Nested approach (“GLUT intensity” ~ Disease_state + (1 | Disease_state:Metadata_SampleID)) to compensate for dependencies within the groups.

3.1 Cytoplasmic area

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: AreaShape_Area ~ Disease_state + (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##       AIC       BIC    logLik  deviance  df.resid 
##  858788.8  858824.0 -429390.4  858780.8     48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.8402 -0.7599 -0.2039  0.6001  3.3928 
## 
## Random effects:
##  Groups                          Name        Variance Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept)  243989   494    
##  Residual                                    2624225  1620    
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                  Estimate Std. Error t value
## (Intercept)        3391.0      221.2   15.33
## Disease_statesPD   -130.4      289.6   -0.45
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: AreaShape_Area
##                Chisq Df Pr(>Chisq)
## Disease_state 0.2028  1     0.6524

3.2 Mean intensity

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: 
## Intensity_MeanIntensity_Corr_GLUT ~ Disease_state + (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##      AIC      BIC   logLik deviance df.resid 
## -54281.8 -54246.6  27144.9 -54289.8    48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -1.5035 -0.7507 -0.4297  0.7451  2.8742 
## 
## Random effects:
##  Groups                          Name        Variance  Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept) 0.0004933 0.02221 
##  Residual                                    0.0191992 0.13856 
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                   Estimate Std. Error t value
## (Intercept)       0.234597   0.009984  23.498
## Disease_statesPD -0.004161   0.013070  -0.318
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: Intensity_MeanIntensity_Corr_GLUT
##                Chisq Df Pr(>Chisq)
## Disease_state 0.1013  1     0.7502

3.3 Intensity cytoplasm ring 10

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: RadialDistribution_FracAtD_Corr_GLUT_10of10 ~ Disease_state +  
##     (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##       AIC       BIC    logLik  deviance  df.resid 
## -116050.2 -116015.0   58029.1 -116058.2     48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -2.9364 -0.7273 -0.1467  0.5683  4.7751 
## 
## Random effects:
##  Groups                          Name        Variance Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept) 0.000109 0.01044 
##  Residual                                    0.005407 0.07353 
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                  Estimate Std. Error t value
## (Intercept)      0.364436   0.004699  77.557
## Disease_statesPD 0.006696   0.006151   1.089
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: RadialDistribution_FracAtD_Corr_GLUT_10of10
##               Chisq Df Pr(>Chisq)
## Disease_state 1.185  1     0.2763

3.4 Intensity cytoplasm ring 9

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: RadialDistribution_FracAtD_Corr_GLUT_9of10 ~ Disease_state +  
##     (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##       AIC       BIC    logLik  deviance  df.resid 
## -158484.7 -158449.5   79246.3 -158492.7     48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -3.8985 -0.7189 -0.0345  0.6786  4.2570 
## 
## Random effects:
##  Groups                          Name        Variance  Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept) 0.0001123 0.01060 
##  Residual                                    0.0022632 0.04757 
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                   Estimate Std. Error t value
## (Intercept)      0.3465913  0.0047517  72.940
## Disease_statesPD 0.0004484  0.0062210   0.072
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: RadialDistribution_FracAtD_Corr_GLUT_9of10
##                Chisq Df Pr(>Chisq)
## Disease_state 0.0052  1     0.9425

3.5 Mean intensity cytoplasm ring 10

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: RadialDistribution_MeanFrac_Corr_GLUT_10of10 ~ Disease_state +  
##     (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##       AIC       BIC    logLik  deviance  df.resid 
## -148962.7 -148927.5   74485.4 -148970.7     48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -5.1728 -0.5275  0.0698  0.6034  4.7720 
## 
## Random effects:
##  Groups                          Name        Variance  Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept) 0.0002095 0.01447 
##  Residual                                    0.0027512 0.05245 
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                    Estimate Std. Error t value
## (Intercept)       0.9570391  0.0064842 147.595
## Disease_statesPD -0.0004615  0.0084894  -0.054
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: RadialDistribution_MeanFrac_Corr_GLUT_10of10
##               Chisq Df Pr(>Chisq)
## Disease_state 0.003  1     0.9566

3.6 Mean intensity cytoplasm ring 9

## Linear mixed model fit by maximum likelihood  ['lmerMod']
## Formula: RadialDistribution_MeanFrac_Corr_GLUT_9of10 ~ Disease_state +  
##     (1 | Disease_state:Metadata_SampleID)
##    Data: data
## 
##       AIC       BIC    logLik  deviance  df.resid 
## -203483.7 -203448.5  101745.9 -203491.7     48736 
## 
## Scaled residuals: 
##     Min      1Q  Median      3Q     Max 
## -5.2625 -0.5655 -0.1075  0.4310  6.3901 
## 
## Random effects:
##  Groups                          Name        Variance  Std.Dev.
##  Disease_state:Metadata_SampleID (Intercept) 1.514e-05 0.003891
##  Residual                                    8.992e-04 0.029987
## Number of obs: 48740, groups:  Disease_state:Metadata_SampleID, 12
## 
## Fixed effects:
##                   Estimate Std. Error t value
## (Intercept)      1.0129478  0.0017539 577.524
## Disease_statesPD 0.0005114  0.0022960   0.223
## 
## Correlation of Fixed Effects:
##             (Intr)
## Diss_sttsPD -0.764
## Analysis of Deviance Table (Type II Wald chisquare tests)
## 
## Response: RadialDistribution_MeanFrac_Corr_GLUT_9of10
##                Chisq Df Pr(>Chisq)
## Disease_state 0.0496  1     0.8238